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ABSTRACT 

We report the discovery of a wide brown dwarf companion to the mildly 
metal-poor ([Fe/H]=-0.24), low galactic latitude (6=1.88°) K4V star HIP 38939. 
The companion was discovered by its common proper motion with the primary 
and its red optical (Pan-STARRSl) and blue infrared (2MASS) colors. It has a 
projected separation of 1630 AU and a near-infrared spectral type of T4.5. As 
such it is one of only three known companions to a main sequence star which 
have early/mid-T spectral types (the others being HN Peg B and e Indi B). Using 
chromospheric activity we estimate an age for the primary of 900±gQQ° Myr. This 
value is also in agreement with the age derived from the star's weak ROSAT de- 
tection. Comparison with evolutionary models for this age range indicates that 
HIP 38939B falls in the mass range 38±20 Mj up with an effective temperature 
range of 1090±60 K. Fitting our spectrum with atmospheric models gives a best 
fitting temperature of 1100 K. We include our object in an analysis of the pop- 
ulation of benchmark T dwarfs and find that while older atmospheric models 
appeared to over-predict the temperature of the coolest objects compared to 
evolutionary models, more recent atmospheric models provide better agreement. 
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1. Introduction 



The identification of a large sample of field L and T dwarfs over the last ~15 years 
has progressed from individual discoveries to large samples of hundreds of objects. Yet i n 



contrast to the tremendous advances made in spectral classification (e.g. iKirkpatrickl 120051 ) . 
accurate determination of the physical properties of ultracool dwarfs remains an unfulfilled 
goal. One inherent difficulty originates from the fundamental nature of substellar objects. In 
the absence of sustained internal energy generation, the temperatures, luminosities, and radii 
for substellar objects of a given mass steadily decline with age — thus diagnosing physical 
properties from photometry and spectroscopy is more complicated than for main-sequence 
stars, which maintain much more constant properties over their lifetimes. In addition, theo- 
retical modeling of emergent fluxes and spectra at such low effective temperatures (<2500 K) 
is challenged by uncertainties in the treatment of dust, clouds, molecules, and chemistry. 

To help overcome these difficulties, the rare subset of brown dwarfs known as "bench- 
marks" is highly valued. This designation is commonly used for any brown dwarf with some 
physical properties much better constrained compared to field objects. Given that three fun- 
damental quantities (e.g., luminosity, temperature, and age) are needed to primarily deter- 
mine a brown dwarf's physical state (ignoring metallicity), a sharper definition wo uld be any 



object for which two fundamental quantities are well-constrained. As described in iLiu et al. 



( 120081 ). substellar benchmarks can be naturally distinguished into "mass benchmarks" and 
"age benchmarks", where the mass or age, respectively, are well-constrained. Mass bench- 
marks are ultracool binaries with directly determined distances (and hence luminosities) and 
dynamical masses. In these binaries, the remaining main properties (temperature and sur- 
face gravity) can be very precisely constrained using evolutionary models. However, given 
the long orbital timescales (>10 yr), the current sample of mass ben chmarks for the L an d 
T spectral classes is limited to five systems (e.g. see compilation in bupuv fc Lhlhoiih l^l 
Ample age benchmarks at hotter temperatures (s pectral types of lat e -M to early-L) have 



been found in you ng (<1-100 Myr) clusters (e.g. iMoraux et al.l l2003t ISlesnick et al.l 12006 



Lodieu et al.l 120071 ) . where the distances and stellar ages are very well-known, and thus the 



1 deacon@mpia.de 

2 There is only one known field system that is both an age and mass benchmark, namely the L4+L4 
binary HD 130948BC, which has a directly dynamical mass determination and a well-constrained age for the 
primary star from gyrochronology ( Dupuv et al. 20091) . 
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luminosites and ages of the brown dwarfs are as well. However, open clusters with ages 
comparable to the field populat ion (>Gyr) are largely dynamically depleted of substellar 



objects (e.g. iBouvier et al.ll2008l ). 



Thus, the search for field benchmarks has largely focused o n brown dwarfs found as 
companions to stars (e.g., see compilation in lFaherty et al.ll2010l ). In this case, the ages of 
the primary stars can also be attributed to their brown dwarf companions, though their typ- 
ical age uncertainties lead to poorer constraints compared to brown dwarfs in short-period 
binaries or open clusters. Three T-dwarf companions have been found at s mall angular sep 



arations through adaptive optics imaging , including the archetype Gl 229B (INakaiima et al. 



1995|), GJ 758B flThalmann et al.ll200a ). and SCR 1845-6357B Oilier et all 120061 ) . Nine 



from wide-field survey data: GL 570D 


(Bureasser et al. 200ol). HD 3651 


3 (Mugrauer et al. 


2006 


). HN PegE 


>(Luhman et al. 


2007 


), Wc 


>lf 940B (Burningham et al. 


2009 


, Ross 458C 


(Goldman et al. 


2010h. HIP 63510C (S 


>cholz 


2010), HIP 73786B 


' Scholz 


2010, 


Murrav et al. 


2011 


). e Indi Bab (Scholz et al. 


2003; 


McCauehrean et al. 2004 


, Gl 337CD 


Wilson et al. 


2001 


Burgasser et al. 2005). Albert et al. (5 


!011) identify CFBDS J022644-062522 as a pos- 



sible wide companion to HD15220 but while the secondary's spectroscopic parallax places it 
at a similar distance to the primary, the proper motion determination on the secondary is 
currently not accurate enough to confirm com panionship. Addition ally, LSPM 1459+0857B 
( IDay- Jones et al.ll201ll ) and WD 0806— 661b ( ILuhman et al.ll201ll ) have been identified as 
wide companions to white dwarfs, the latter of which is much cooler than any known T dwarfs 
but has not yet been spectrally classified. 

Wide-field surveys represent fertile ground for identifying benchmarks. Most wide 
benchmark companions have resulted from large-area searches for free-floating ultracool 
dwarfs, which are then realized to be comoving with higher mass stars. More limited efforts 
have been car ried out to expressly identify wide benchm arks around well-defined samples of 
primary stars (jPinfield et al.ll2006l ; IDay- Jones et al.ll2008l ). With a sufficiently large sample of 
benchmarks, it may be possible to provide a direct empirical calibration that ties the under- 
lying physical parameters (e.g., age and temperature) to observable brown dw arf properties 



e.g. IPinfield et al.ll2006f ). To this end, the Pan-STARRSl (PS1) Telescope flKaiser et al. 



20101 ) is a welcome addition to help boost the census of benchmark brown dwarfs. 



Situated on the summit of Haleakala on the island of Maui in the Hawaiian Islands, 
PS1 will be the leading multi-epoch optical survey facility over the next several years. PS1 
is the first of four planned 1.8- meter tel escopes that will c omprise the Panoramic Survey 
Telescope And Rapid Response System (IKaiser et al.l 120021 ) . Full science operations with 
PS1 began in May 2010, and fabrication of the second telescope (PS2) has already begun. 
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So far, c ommissioning and regular surv ey data have been used to search for Trans-Neptunian 
Objects (W ang et al.ll2010h . T dwarfs ( IDeacon et al.ll201ll . lLiu et al.ll2011al ) and supernovae 



(e.g. 



Botticella et all l2010h 



The telescope is executing several astronomical surveys of 
varying depth, area, cadence, and filter complements. The most interesting of these for 
discovering free-floating and companion brown dwarfs is the 2>tt Survey, which is scanning 
the entire sky north of S = —30 ° (3tt steradians) in five filters (gpi, rp±, ip±, zp\ and ypi, 
Stubbs et al.ll2010l . iTonryi l201ll in prep.) at six separate epochs over 3.5 years, with each 
epoch consisting of pairs of exposures taken f»25 min apart. These multiple epochs and the 
use of the 0.95-1.03 /im ypi band make the survey well-suited for ide ntifying nearby cool 



objects such as brown dwar fs by their proper motions and parallaxes ( IMagnier et al 
Beaumont fc Magnierll2010l ). 



2008 



In this paper, we present the discovery of a T4.5 dwarf around the K4V star HIP 38939 
found from a dedicated search for wide late-L and T dwarf companions using the PS1 and 
2MASS. In § 2, we describe our mining of these datasets to search for companions to main- 
sequence stars with proper motions from Hipparcos. Spectroscopic confirmation is presented 
in § 3. In § 4 we examine the physical properties of our new discovery, along with a look at 
the total compilation of benchmark T dwarfs now known. 



2. Identification in PanSTARRSl + 2MASS data 



Pan-STARRSl began full survey operations in May 2010. Hence for a given filter and 
location on the sky, the current 3tt survey data consist of only one pair of images taken at 
the same epoch or multiple pairs separated by a short time baseline. So we can use these 
data for studies of nearby brown dwarfs, we have undertaken a p roper motion survey by 
combining PS1 data with t hose from 2MASS (pkrutskie et al.ll2006f ). This method was used 
to identify new T dwarfs in IDeacon et al.l ( 1201 ll ). We constructed a database of high proper 
motion objects containing all PS1 37r survey commissioning and survey data taken from the 
beginning of February 2010 to the end of January 2011. This was done by identifying PS1 
ypi objects with no 2MASS detection within one arcsecond. We required that the PS1 object 
was detected on at least two images, was classified as a good quality point source detection, 
and had a photometric error less than 0.2 magnitudes in the yp\ band. The double detection 
requirement excludes contamination due to asteroids and reduces the number of instrumental 
artefacts. These objects are then paired with 2MASS J-band point sources which had an A, 
B or C photometric quality flag designation (corresponding to S/N>5), which were within 
28" and which had no corresponding PS1 detection. This pairing radius was chosen to allow 
the detection of objects with proper motions up to 2"/yr given the 14-year maximum epoch 



-5 - 



difference between PS1 and 2MASS. 



To identify companion objects to nearby stars, we ran queries on ou r proper motion 
database around all stars in the Hipparcos catalogue (jPerryman et al.lll997| ). We chose stars 
with proper motions above 0.1"/yr and parallaxes more significant than 5a. We searched for 
companions with projected separations less than 10,000 AU at the distance of the Hipparcos 
star. Each candidate was further filtered to have proper motion vectors within 5cr of the 
primary assuming a nominal accura cy of 16 milliarcsecond s per yeai§| in each coordinate for 
our proper motions and errors from iPerryman et al.l ( 119971 ) for t he primary W e then picked 
out candidate ultracool dwarfs using the color selections found in lDeacon et al.l (1201 if ) which 
are sensitive to late-L and T dwarfs. These are ypi — J > 2.2, J — H < 1.0 and zp^ — J7pi > 
0.6 with no c orresponding R or / band detections in SuperCOSMOS ( IHambly et al.l l200lh 
or USNO-B ( IMonet et al.ll2003l ) brighter than 20.5 and 18.5 magnitudes respectively. No 
Galactic latitude cut was applied. Taking into account area lost due to gaps between chips in 
the PS1 camera and areas of multiple coverage where these gaps can be filled in, we searched 
16,900 sq.deg. After examining 1008 candidates by eye to remove spurious associations and 
objects with clear but uncatalogued companions in the USNO-B photographic plate images 
380 candidates remained. Of these, our clearest late-type candidate was a companion to 
HIP 38939 (a.k.a. CD-25 5342, GJ 9246, HD 65486, NLTT 18729, SAO 174889). We have 
not yet completely followed up our remaining candidates, hence this work to date should not 
be considered a complete survey for late-type companions to Hipparcos stars. 

Details of the companion are shown in Table [TJ Its proper motion is a good match 
to that of the prim ary, differing by 1.5er. At the distance to the primary (18.5±0.4 pc; 
van Leeuwenl 120071 ) the angular separation of 88 arcseconds means a projected separation 
of 1630 AU. Finder charts are shown in Figure [TJ The relatively dense field is due to the 
low Galactic latitude of HIP 38939 (b = 1.88°). This explains why the object has gone 
undiscovered until now, as most brown dwarf surveys have avoided the Galactic plane. 



3 This is a conservative estimate based on a comparison between the measured proer motions for stars in 
the lLepine fc Sharal |2005l ) catalog and the proper motions we derived for these objects. 
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9pi r pi 'pi z pi Ypi 




J H K s W1 W2 

Fig. 1.— Images of HIP 38939B from Pan-STARRSl (g P1 , r P1 , i P1 , z P1 , y P1 ), 2MASS (J, 
H, K s ) and WISE (Wl, W2). The white circles are centered on the position in the Pan- 
STARRSl y p i image (solid line) and the 2MASS position (dashed line). Images are one 
arcminute across with North up and East left. The cutout images in the upper left of the 
Pan-STARRSl panels are 5" across and centered on the position in the Pan-STARRSl y P1 
image shown. 
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Table 1. The HIP 38939 system. 



HIP 38939A HIP 38939B 



RA (J2000) 


07 58 04.13 a 


07 58 01.61 1 


Dec (J2000) 


-25 37 33.7 a 


-25 39 01.4 y - 


H a co&{5) ("/yr) 


0.362±0.001 a 


0.353±0.006 


^ (7yr) 


-0.245±0.001 a 


-0.244±0.006 


<) 


0.054±0.001 a 




Spectral Type 


K4V C 


T4.5V d 


B (mag) 


9.50° 




V (mag) 


8.44 c 




zpi (AB mag) 




20.17±0.15 f 


y P1 (AB mag) 




18.66±0.08 f 


2MASS J (mag) 


6.51±0.02s 


16.12±0.08 s 


2MASS H (mag) 


5.94±0.02s 


15.80±0.12s 


2MASS K s (mag) 


5.83±0.02s 


>15.86 g 


M bo i (mag) 


7.93 h 


17.04±0.18 d 


r e // (K) 


4683±35 h 


1090±gf 
1100 j 


logi # (cm/s 2 ) 


4.5 h 


5.0=1=0.3* 
4.5 1 



Synthesized colors^ 



J - H (2MASS) (mags) 0.139±0.003 

H - fsT s (2MASS) (mags) -0.104±0.006 

J - fT s (2MASS) (mags) 0.03±0.006 

Jmko ~ J2MASS (mags) -0.217±0.007 

J - if(MKO) (mags) -0.127±0.003 

H — K(MKO) (mags) -0.191±0.007 

J - KiMKO) (mags) -0.318±0.006 
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All astrometry f or the primary comes from 



van Leeuwenl ( 120071 ) 



b Epoch 2010.0 



Gray et all (120061 ) 



d This work 



Koen et all (12010 ) 



Pan-STARRSl Image Processing Pipeline. 
Zeropoints are p rovisionally calibrated relative 



to the 2MASS (ISkrutskie et al J 12009). Tvcho 
jjfog et al.ll2000h and USNO-B (Monet et al 



20031 ) surveys. 



g Skrutskie et al. 



(120061 ) > indicates the 
2MASS magnitude is a 97% confidence upper 
limit. 



ICasagrande et all ( EOlOn 

'This work, evolutionary models, see Section 



JThis work, atmospheric models, see Section 
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2.1. Data from other surveys 



HIP 38939B does not appear in Data Release 8 of the Sloan Digital Sky Surve y (Aihara et al. 



2011 ) and is too f ar south to appea r in the UKIDSS Galactic Plane Su rvey (ILucas et al. 



20081 ). the IPHAS flDrew et al.ll2005h or UV EX surveys flGroot et al.ll2009f ). The object falls 
in the area of the WISE (I Wright et al.ll2010l ) Preliminary Data Release but has no detection 
in the catalogue. However visual inspection indicated there was a bright H / 2-band object 
at the PS1 posi tion of the source. We r etrieved FITS images from the WISE database and 
ran SExtractor (IBertin fc Arnoutslll996l ) on the area. We then used the WISE catalog mag- 
nitudes for other objects in the field to calculate the zero points of the images. The W2 
magnitude for this HIP 38939B is 13.82 ± 0.04 mag, with a distance of 0.6" from the 2010.0 
position. In the Wl-band there is a faint (Wl = 15.54 ± 0.13 mag) object 4.8" away from 
the PS1 position which appears to be blended with another object of similar magnitude (see 
Figured]). Hence we could not determine an accurate Wl-band magnitude for HIP 38939B. 



2.2. Likelihood of chance alignment 

We searched for companions around 9363 nearby Hipparcos with projected separations 
less than 10,000 AU. The sum of these search areas around our stars comes to 260 sq.deg. 
The dwarfarchives.org website reports 90 T dwarfs known to an approximate 2MASS depth 
of J ~ 16.5 (about 1 per 400 sq.deg.). Hence there is a reasonable chance that one field 
T dwarf will appear within our search area by coincidence. The probability of a random 
object having the same proper motion as HIP 38939A was estimated by taking objects 
from our high-speed Pan-STARRS1-2MASS catalogue within one degree of the primary and 
estimating what fraction had a similar proper motion to the primary. Only 5 out of 31349 
objects had proper motions within 5a. Hence we believe (assuming the velocity distribution 
of T dwarfs is not substantially different from that of objects in the HIP 38939 field) that 
the probability of an unassociated T dwarf having the same proper motion as the primary 
is <0.02%. 



3. Spectroscopy 



Spectroscopic c onfirmation of HIP 38939B was obtained using the SpeX spectrograph 
( iRayner et al.l 120031 ) on the NASA Infrared Telescope Facility on Mauna Kea on 31 March 



4 http: / /dwarfarchives.org 
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2011 UT. Observing conditions were clear with good seeing. We used the low resolution, 
single-order prism mode to obtain a 0.8 — 2.5 fim spectrum. Observations were made at an 
airmass of 1.89 using the 0.5x15 arcsecond slit, yielding a spectral resolution of ~120. The 
total integration time was 960 seconds. The slit was oriented at the parallactic angle to min- 
imize atmospheric dispersion. We nodded the telescope in an ABBA pattern with individual 
exposures of 120s. The spectra were extracted, wavelength calibrated a nd telluric corrected 



using a contempo raneously observed A0V star in the SpeXtool package (jCushing et al.ll2004 
Vacca et al.l 120031 ) . The final reduced spectrum is in Figure [2J 



HIP 38939B 




1.2 1.4 1.6 1.8 2.0 

Wavelength (microns) 



2.4 



Fig. 2. — The near- infrared spectrum of HIP 38939B. The absorption bands characteristic 
of a T dwarf can clearly be seen. Note the spectrum has been normalized to its maximum 
value. 
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3.1. Spectral Classification 



We assigned a near-IR spectral type for HIP 38939B us ing the flux indices of lBurgasser et al. 



(l2006bl ) and the polynomial relations of iBurgasserl (j2007al ). The spectrum was also visually 



i nspect ed by comparing with IRTF/SpeX prism spectral standards from iBurgasser et al. 



(l2006bl ). In this process the spectrum was normalized to the peak fluxes of the standards in 
the J, H and K bands individually and the depth of the H2O and CH4 absorption bands 
were examined. These spectral indices are shown in Table [2j From visual comparison, we 
find the best fitting spectral type is T5: our spectral index measurements indicate a type of 
T4.5. We allow our spectral index measurement to take precedent over the visual comparison 
and we assign a spectral type of T4.5. 

From our spectrum we synthesized a color term be tween the MKO and 2 MASS J bands. 
This gave a J band magnitude in the MKO system (ITokunaga et al.l 120021 ) of 15.90±0.08 
mag. We then used this combi ned with an aver age of the bright and faint spectral type to 
absolute magnitude relations in iLiu et al.l ( 120061 ) to estimate a photometric distance to HIP 
38939B of 20.8^21 P c - This differs from the Hipparcos dist ance to the primary of 18.5±0.4 
pc by less than la. Using the similar distance estimator of iGoldman et al.l (120101 ) we get a 
distance of 20.8±4.3 pc, within la of the distance to the primary. 



4. Discussion 



4.1. The characteristics of HIP 38939 A 



Gray et al.l (120061 ) classify HI P 38939A as a K4V. Th is is consistent with the value of 



T, 



eff 



Casagrande et al.l (120101) . who also estimated a sub-solar 



Santos et al.l J2OO5I ) measured 4660 ± 66 K and [Fe/H] 



4683 ± 35 K measured by 
metallicity of [Fe/H]=-0.24, while 
0.33 ± 0.07. The Hipparcos satellite measured a trigonometric parallax of 54.91 ± 1.15 mas 
and proper motions of fi a cosS = 362.78 ± 0.070 mas/yr and fj,s = —245.89 ± 0.63 mas/yr 
( Ivan Leeuwenll2007l ). 



In order to properly characterize HIP 38939B we must constrain its age by estimating the 
Table 2. Spectral indices for HIP 38939B. 



H20-J 


CH4-J 


H2O-H 


CH4-H 


CH4-K avg/RMS Visual Final 


0.303 (T4.5) 


0.492 (T3.6) 


0.388 (T4.4) 


0.453 (T4.9) 


0.244 (T4.8) T4.4±0.5 T5 T4.5±0.5 
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age of its prima ry. One of the key estimators for the ages of F, G and K type star is chromo- 
spher ic activity ( Wilson 1963 . Soderblom et alll991 . Donahue 1998 . Mamajek fc Hillenbrand 
20081 ). This is the by-product of the faster rotation rates of younger stars which drive cor- 
respondingly more active dynamos and hence higher chromospheric activity. The chromo- 
spheric activity is characterised by the emission in the calcium H and K lines and is typically 
expressed as, 

R'hK = R-HK — Rphot (1) 

Where R p hot is the photospheric contribution to the flux in the Calcium H and K lines, and, 



R 



UK 



1.34 x 10~ 4 C cf S 



(2) 



Here C c f is a correction factor dependent on color (IMiddelkoopI 1 1 9821 . iNoyes et al.lll984f ) and 
S is a ratio of spect ral indices in the H and K lines to continuum indices as defined by 
Vaughan et "all f)l978h . 



Gray et al.l ( 120061 ) classify HIP 38939 A as bei ng "active" meaning it has a log 
value between —4.75 and —4.2. IWright et al.l (120041 ) measure the mean S value of the object 
to be 0.757, but due to concerns about calcu l ating the correction factor do not measure 
an age. We use t he relations in iNoves et al.l (119841) t o cal culate that C c f = 0.314 based 
on B — V — 1.06 ( Koen et al. 201oh . As~ Noyes et al. ( 1984 ) suggest that the photospheric 
contribution to the Calcium H and K emission for stars redder than B — V = 1.0 is negligible, 
we assume R p hot = 0. Hence we derive a log 10 R' HK value of —4.5. 

19981) 



Co mparing to indivi dual cluster datasets from the 680 Myr old dPerrvman et al. 



Hyades ( jSo derbloml 1 1 98 5l ) and the 500 Myr old Ursa Major moving group (IKing et al. 



20031 ). 



we find that HIP 38939A has chromospheric emission consistent with the ages of these 
two young clusters/associations. However we find that our activity value is no t a good 
fit when compared to dat a from the younger (100-120 Myr iMartm et al.l Il998l ) Pleiades 
dSoderblom fc Mavor|[l9~93T ). 



The most current age-activity-rotation relation is that of lMamajek fc Hillenbrand! (120081 ) 
which unf ortunately does not apply to stars redder th an B — V = 0.9 . Using the earlier 



relation of iDonahud ( 119981 ) we derive an age of 870 Myr. iDonahud (119981 ) estimates a spread 
of 0.5 dex for the age of the Sun from a single measurement of the chromospheric activity at 
any point across the solar cycle. This value is consistent with the re sults for younger, co eval 
binaries and is larger than the intrinsic variability of 13% quoted by lWright et al.l (120041 ) for 



their S values. Hence we assign errors to the Donahue value of 900±Jgg° Myr. 

Another indicator of stellar a ctivity and hence youth is X-ray emission. We searched the 
ROSAT Faint Source Catalogue ( IVoges et al.ll2000l ) and identified an extremely faint source 
18" away from HIP 38939A. This distance is slightly larger than the quoted ROSAT la 
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positional uncertainty for this source of 15". As there is always the possibility of association 
with an unrelated source, we carried out the following calculation. There are 25 ROSAT 
sources within 3 degrees of the primary. This corresponds to a surface density of 0.88/sq.deg. 
If we assume a ROSAT source has to be within 30" of the primary to be associated then the 
probability of a chance alignment is 0.02%. Hence we assume the X-ray source is associated 
with HIP 38939A. The source has X-ray counts of (3. 66 ± 1.48) x 1 ~ 2 s~ l and a hardness 
ratio of HR1 = —0.76 ± 0.27. Using the relations of Schmitt et al.l (1995 ) we cal culate an 
X-ray luminosity of L x = (6.43 ± 3.37) x 10 27 ergs/s. ICasagrande et al.l (120061 ) derive a 
bolometric luminosity of 0.175 L & for HIP 38939A. From this we calculate a fractional X-ray 
luminosity R x = Ior(L x /L hn] ) = log(Y9.55 ± 5.00) x 10~ 6 ). Using the X-ray luminosity - age 
relation from lMamajek fc Hillenbrand! (120081 ). we derive an age of 900±gg{j° Myr taking in to 
account both the measurement errors on the X-ray flux and the quoted 0.4 dex uncertainty 
in the age relation. While this is a low significance X-ray detection, we note that the X-ray 
derived age is in agreement with those derived from calcium H and K emission. 



4.2. Comparison of HIP 38939B with evolutionary models 



We derived the bolometric luminosity for HIP 38939B by converting the 2MASS J- 
band magnitude into the MKO system using a color term synthesized from our spectrum, 
converting to absolute magnitude using the Hipparcos parall ax for the primary and applying 
the bolometric correction relation found in iLiu et al.l (J2010J). This resulted in a bolometric 
magnitude of M\ M \ = 17.04 ± 0.18 mag. We then compared the range of ages derived for th e 
primary along and our calculated M bol to the evolutionary models of iBurrows et al.l (119971 ). 
Figure [3] shows where the secondary lies in comparison with these evolutionary models. 
Table [3] also shows the values of various parameters for different ages. Additionally we ran 
a Monte Carlo simulation with our calculated bolometric magnitude and with a uniformly 
distributed range of ages derived from the primary. The results for these are also shown in 
Table [3J 



4.3. Comparison with atmosphere models 



Atmospheric models provide an independent means of derivi ng the physical pro perties 
of HIP 38989B. We fit the solar metallicity BT-Settl-2010 models (Ullard et al.lboioh to our 



SpeX spectrum following the procedures in ICushing et al.l (120081 ) and iBowler et al 



(120091 ). 



The grid of models spans effective temperatures between 500-1500 K (AT c g=100 K) and 
gravities between 10 4 '°-10 5 ' 5 (cgs; Alog g=0.5). The SpeX spectrum was first flux calibrated 



500 1000 1500 2000 

T e „ (K) 



Fig. 3. — Our derived age and bolom e tric l uminosity for HIP 38939B plotted against the 
evolutionary models of IBurrows et al.l (Il997l ). The blue-colored hatched areas indicate the 
constraints on the parameters given the uncertainties on Mbd and the age of the primary. 
The solid lines represent isochrones with the annotated values of log(age). The two unmarked 
isochrones at the bottom are log(age)=8.0 and 7.7. The dashed lines are lines of equal mass 
with the appropriate values noted on the right-hand side of the plot. See Table [3] for more 
details on the individual masses for each particular age and the Monte Carlo simulated 
parameters. 
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to the 2MASS J-band photometry and then fit between 0.8-2.4 /mi in a Monte Carlo fashion. 
The 1.60-1.65 /mi region was excluded from the fits because the methane line list in the 
models is incomplete at those wavelengths (e.g., Saumon et al. 2007). For each Monte Carlo 
trial we alter the flux calibration scaling factor and SpeX spectrum by drawing random values 
from normal distributions based on photometric and spectral measurement uncertainties, 
respectively. We then fit the grid of models to each artificial spectrum and repeat the 
process 10 3 times. 

The best-fitting BT-Settl-2010 model has T eff =1100 K and log g=4M This is in good 
agreement with the evolutionary model predictions of T o fj=1090lgQ K and log g=5.0±0.3 
especially considering the coarse gridding of the atmospheric models. In Figure H] we show 
the flux-calibrated SpeX spectrum along with photometry from PS1 {yp\ and Zp\ bands), 
2MASS (J, H, and an upper limit on K$), and WISE (W2). The monochromatic flux 
densities agree well with the spectrum between 0.8-5 /im. The best-fitting atmospheric 
model provides a good match to the spectrum except at ~1.15 /mi, 1.65 /im, and 2.15/an, 
which correspond to absorption features from CH4. 

Using our SpeX spectrum together with the J-band photometry measured by 2MASS, 
we also synthesize photometry of HIP 38989B in the 2MASS H and K s bands and our 
spectral measurement errors. The synthesized magnitudes for 10 4 Monte Carlo trials are 
/J=15.98±0.08, and i^5=16.09±0.08 mag. The H band synthetic magnitude is within 2a of 
the 2MASS measurement, which is near the detection limit of the survey. The synthesized 
colors are J - JJ=0.139±0.003, H - AT S =-0.104±0.006 and J - if g =0.03±0.006 mag. The 
uncertainties in the synthesized colors are much smaller than the magnitudes because they 
only include spectral measurement uncertainties. 



4.4. The population of benchmark brown dwarfs 

As noted in Section 1, brown dwarf binaries can be used to test atmospheric and evo- 
lutionary models of substellar objects. Atmospheric models of brown dwarfs will not fully 
match the shape of observed spectra due to missing opacity lists for H 2 0, CH 4 and NH 3 



(ILeggett et al.l 120071 ). In order to produce a more complete picture of how these models 
compare with the growing population of benchmark brown dwarfs, we compiled a list of 
all T dwarf companions with effective temperatures derived from both model atmospheres 



5 We also carried out a fit to the models including the 1.6-1.65 [im spectral region. The best fit model 
was T e ff=1200 K/log g=5.5, which differs slightly from the 1100 K/4.5 dex best fit when that wavelength 
region is excluded. 
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Fig. 4. — Left: The comparison between the observed SpeX spectrum (black) and photome- 
try (green circles for data points, green triangles fo r upper limits) and the best fitting model 
atmosphere (red) from the BT-Settl-2010 models (jAUard et al.ll2010l ). The spectrum is flux 
calibrated using the J-band photometry from 2MAS S. The magnitude to flu x conversions 



used the zeropoints and effective wavelengths from iRieke et al.1 (120081 ) and IWright et al. 



(120101 ) for 2MASS and WISE respectively, and the PS1 image processing pipeline AB zero- 
points and the effective wavelengths from Tonry (in prep), for the PS1 data. The horizontal 
error bars represent the filter bandpass. The WISE 4.6 /im point is consistent with the 
best-fit model atmosphere, which assumes chemical equilibrium (no vertical mixing). Upper 
right: The xl contours for the comparison between the models and the observed spectrum, 
the contours represent values of 70, 100, 150, and 250. The best fitting model (1100 K, 
log (7=4.5) is shown as a yellow star with the blue star and error bars representing the r esults 
of Monte Carlo calculations based on the evolutionary models of iBurrows et al.l (119971 ). The 
two agree within the stated error in T e ff and differ by only one grid-point in logg. When 
the region of the spectrum from 1.6-1.65 /im (previously excluded due to incomplete model 
line lists) is included, the best fitting gravity moves to log g—5.5. This is still within one grid 
point of the evolutionary model calculations. Lower right: T he different radii as a fraction of 
the radius of Jupiter from the model atmosphere fitting.. See lBowler et al.l (120091 ) for details 
on this method. 
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and evolutionary models. In order to be included, an object had to have an effective tem- 
perature derived from its bolometric luminosity and a separate temperature derived from 
model atmosphere fits to the observed spectrum. These are listed in Table H] along with 
the models used for each calculation. The two temperatures for each object are plotted 
against each other in Figure [51 Evolutionary models depend on atmospheric models to set 
their boundary conditions. In many of the comparisons we make, the evolutionary model 
boundary atmosphere is different from that used to derive the atmospheric effective tem- 
perature. In cases where the same atmospheric model is used this has been noted in both 
Figure [5] and Table HI Many of our objects have evolutionary model effective temperatures 
derived from different sou rces. Table H] shows th at t here is no signi f icant diffe rence between 



the de rived temperatures. iLuhman et all (120071 ) use Burrows et all ( 119971 ) and iBaraffe et al. 



( 120031 ) evolutionary models to derive effective temperatures for HN Peg B and HD 3651 
B finding that there is no significant difference between the two. There appears to be a 
bias towards lower temperature objects having atmospheric temperatures which are higher 
than those derived from evolutionary models. It should be note d that the best fits for these 
objec ts are those from the mos t recent models, BT-Settl 2 010 (Allard et al.ll2010l) for Ross 
458C iBurningham et al l ( boilh . fits to lSaumon fc Marlevi (fioosh bvlBurgasser et all (hoioh 
and a set of model s calculated by Saumon for Leggett et all ( 2010 Vs study of Wolf 940B. 



Dupuy et all ( 120101 ) noted that atmosphere models for late-M dwarfs also predict higher 



temperatures than the evolutionary models. 
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Table 3. Derived properties for HIP 38939B using our inferred age range for the primary. 



logio(age) 
yr 


Mass 
M Jup 


T eff 

K 


log(#) 
cm/s 2 


8.40 


19 


1020 


4.6 


8.70 


27 


1060. 


4.8 


8.90 


34 


1090 


5.0 


8.95 


36 


1090 


5.0 


9.00 


38 


1100 


5.0 


9.05 


40 


1110 


5.1 


9.10 


42 


1120 


5.1 


9.30 


52 


1150 


5.3 


9.40 


59 


1170 


5.3 



Monte Carlo Results 



8.90±0.5 34+}? 1090+12 5.0±0.3 



Table 4. Known T dwarf companions to main sequence stars with effective temperatures from model atmospheres 

(T e ff,atm) and evolutionary models (T effievo i). 



Object 


SpT 
K 




Evolutionary Model 




Reference 
K 


-^"e/ f,atm 


Atmosphere Model 


Reference 


HIP 38939B 


T4.5 a 


1090±™ 


Burrows et al. (1997) 




This work 


1100 


Allard et al. (2010) 


this work 


e Indi Ba 


Tl b 


1368±17 c 


Baraffe et al. (2003) 




King et al. (2010) 


1320±20 


Allard et al. (2010) 


King et al. (2010) 














1275±25 


Burrows et al. (2006) 


KasDer et al. (2009) 


e Indi Bb 


Tg b 


993±18 c 


Baraffe et al. (2003) 




King et al. (2010) 


910±30 


Allard et al. (2010) 


King et al. (2010) 


HN Peg B 


T2.5 d 


1065±50 


Leggett et al. (2008) 




Leggett et al. (2008) 


900±25 
1115 


Burrows et al. (20061 
Leggett et al. (20081 


Kasper et al. (2009) 
Leggett et al. (2008)* 






1130±70 


Burrows et al. (1997) + Baraffe et al. 


(2003) 


Luhman et al. (2007) 






Gl 570 D 


T7.5 C 


794±60 


Burrows et al. (19971 




Geballe et al. (2001) 


800±20 


Burrows et al. (2006) 


Burgasser et al. (2006a) 














900±50 


All l j_ l /nriAi \ 

Allard ct al. (2001) 


Liu et al. (20llb) 














800±50 


Allard & Frevtag (2010) 


Liu et al. (201lb) 


HD 3651 B 


T7.5 d 


810±30 
810±50±70 


Burrows et al. ( 1997) 

Burrows et al. (1997) + Baraffe ct al. 


(2003) 


Liu et al. (2007) 
Luhman et al. (2007) 


810±30 
800±50 


Burrows et al. (2006) 
Allard et al. (2001) 


Burgasser (2007b1 
Liu et al. (2C)llb)* 


Ross 458 C 


T8.5p f 


695±60 


Saumon & Marlev (2008) 




Burniiigliam ct al. (201T) 


850±50 
725±25 


Allard & Frevtag (2010) 
Allard et al. (2010) 


Liu et al. (2011b) 
Burninghanf^t al. (2011) 






650±25 


Saumon & Marlev (20081 




Burgasser et al. (2010) 


635 ±30 


Saumon & Marlev (2008) 


Burgasser et'al. (2010) h * 














760±J§ 


Saumon & Marlev (2008) 


Burgasser et al. (2010V * 














900±50 


Allard et al. (2001) 


Liu et al. (2011b) 














850±50 


Allard & Frevtag (2010) 


Liu et al. (2011b) 


Wolf 940 B 


T8.5 1 


590±40 


Saumon fc Marlev (2008) 




Leggett et al. (2010) 


605±20 


Leggett et al. (2010) 


Leggett et al. (2010)* 






570±25 


Baraffe et al. (2003) 




Burniiigliam ct al. (2009) 


800±50 


Allard et al. (2001) 


Liu et al. (2011b)* 














650±50 


Allard & Frevtag (2010) 


Liu et al. (2011b) 





















* Evolutionary model boundary atmosphere model is the same as the model used to derived the atmospheric effective temperature. 
a This work 

,: McCaughrcan ct al . ( 2001) 

c This is based on an age of 4.0±0.3. A number of different indicators give a wide range of ages for the e Indi system. An older age would bring the evolutionary models into better 
agreement with the atmosphere models llKing et al.ll20ld . lLiu et al.ll2010l 

jLuhman et alj J2007h 

^Burgasser et al. (2006b) 



'BurniiiKham ct al. 
g cloudy model 
h cloudless model 
'Burningham et al. (2009) 



to 
O 
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5. Conclusions 

We have discovered a T4.5 companion to the nearby K4V star HIP 38939 using Pan- 
STARRS1 and 2MASS data. We estimate the age of the primary as 900±gQQ° Myr from 
which we deduce that the companion has a mass in the range 38±20 Mj up . The effective 
temperatures and surface gravities derived from model atmosphere fitting and evolutionary 
models are consistent. Given its close proximity to Earth (18.4 pc), this object will be a 
high priority target to identify if it itself is a binary. Such obi ects are of great importance 



for testing physical models of brown dwarfs (ILiu et al.l 120081 ) . This discovery indicates Pan 



STARRS 1 data can be used to explore crowded regions of the sky such as the Galactic 
plane often ignored by other studies. In the future the wide-field, multi-epoch nature of the 
Pan-STARRSl 3n survey will lead to a well-defined sample of wide substellar companions 
to nearby main sequence stars. 



The PS1 Surveys have been made possible through contributions of the Institute for 
Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck 
Society and its participating institutes, the Max Planck Institute for Astronomy, Heidel- 
berg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hop- 
kins University, the University of Durham, the University of Edinburgh, Queen's University 
Belfast, the Harvard-Smithsonian Center for Astrophysics, and the Los Cumbres Observa- 
tory Global Telescope Network, Incorporated, the National Central University of Taiwan, 
and the National Aeronautics and Space Administration under Grant No. NNX08AR22G is- 
sued through the Planetary Science Division of the NASA Science Mission Directorate. They 
would like to thank Dave Griep for assisting with the IRTF observations. This research 
has benefitted from the SpeX Prism Spectral Libraries, maintained by Adam Burgasser at 
http:/ /www.browndwarfs.org/spexprism, This publication makes use of data products from 
the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts 
and the Infrared Processing and Analysis Center /California Institute of Technology, funded 
by the National Aeronautics and Space Administration and the National Science Founda- 
tion. This research has benefitted from the M, L, and T dwarf compendium housed at Dwar- 
fArchives.org and maintained by Chris Gelino, Davy Kirkpatrick, and Adam Burgasser. This 
research was supported by NSF grants AST-0507833 and AST09-09222 (awarded to MCL), 
AST-0709460 (awarded to EAM), AFRL Cooperative Agreement FA945 1-06-2-0338, and 
DFG- Sonderforschungsbereich 881 The Milky Way. Finally, the authors wish to recognize 
the very significant cultural role that the summit of Mauna Kea has always had within the 
indigenous Hawaiian community. We are most fortunate to conduct observations from this 
mountain. 

Facilities: IRTF (SpeX), PS1 
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Fig. 5. — The comparison of the evolutionary and atmospheric model effective temperatures 
for the objects listed in Table HI The error bars represent the range of values quoted in 
the references and the dashed line marks a one to one agreement between the two tempera- 
tures. The gray bars link the datapoints of objects with multiple values of the atmospheric 
effective temperature. HIP 38939B is marked by a s tar. The "othe r models" plot t ed ar e 



model calcula t ions b ased on ISaumon &: Marlevl (120081 ) carried out for iLeggett et aU ( 120101 ) 



Leggett et al.l (120081 ) and iBurgasser et al.l (l2010l ). Square symbols represent objects where 



the evolutionary model effective temperature is derived from a model which uses the same 
atmospheric model as a boundary condidtion as was used to derive the atmospheric model 
effective temperature. It appears that for the coolest T dwarfs, only the most recent models 
provide good agreement between the two temperatures. For the components of e Indi B the 
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